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Abstract 

We describe algorithms which address two classical problems in lattice geometry: the lattice cov- 
ering and the simultaneous lattice packing-covering problem. Theoretically our algorithms solve the 
two problems in any fixed dimension d in the sense that they approximate optimal covering lattices and 
optimal packing-covering lattices within any desired accuracy. Both algorithms involve semidefinite 
programming and are based on Voronoi's reduction theory for positive definite quadratic forms, which 
describes all possible Delone triangulations of Z d . 

In practice, our implementations reproduce known results in dimensions d < 5 and in particular solve 
the two problems in these dimensions. For d = 6 our computations produce new best known covering as 
well as packing-covering lattices, which are closely related to the lattice Eg. For d = 7, 8 our approach 
leads to new best known covering lattices. Although we use numerical methods, we made some effort 
to transform numerical evidences into rigorous proofs. We provide rigorous error bounds and prove that 
some of the new lattices are locally optimal. 
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1 Overview 



Two classical problems in the geometry of numbers are the determination of the most economical lattice 
sphere packings and coverings of the Euclidean ci-space K d . In this paper we describe algorithms for the lat- 
tice covering and the simultaneous lattice packing and covering problem (lattice packing-covering problem 
in the sequel). 

Roughly speaking, both problems are concerned with the most economical way to cover E d . In the case 
of the lattice covering problem, the goal is to maximize the volume of a fundamental domain in a lattice 
covering with unit spheres. Roughly speaking, we want to minimize the number of unit spheres which 
are needed to cover arbitrarily large but finite regions of K d . The objective of the lattice packing-covering 
problem is to maximize the minimal distance between lattice points in a lattice covering with unit spheres. 

The aim of this paper is to give an introduction to the mathematical tools that allow us, at least in theory, 
to solve the two problems computationally. For a fixed dimension d, our algorithms approximate optimal 
covering lattices and optimal packing-covering lattices within any desired accuracy. In this overview we 
want to describe the structure of the paper. 

The basic concepts and notations, which we use throughout the paper, are in Section 13 There we also 
give a precise definition of the two problems under consideration. The reader familiar with sphere packings 
and coverings, as well as lattices and their relation to positive definite quadratic forms, may skip this section. 

In Section |3] we review known results and the history of the lattice covering problem, and in Section |4] 
we review known results and the history of the packing-covering problem. 

Our algorithms as well as the known results by other authors are mainly based on a reduction theory 
for positive definite quadratic forms by Voronoi. We give a detailed description of this main ingredient in 
Section|5]with a special focus on computational implementability. 

The other main tool comes from convex optimization theory. Semidefinite programming problems and 
determinant maximization problems are briefly described in Section |6] We describe how duality theory 
together with rational approximations can be used to provide rigorous error bounds. Both problems have in 
common that one has to minimize a convex function on variables that satisfy some linear matrix inequalities 
(LMIs). 

In Section we describe how the constraint that a lattice gives a unit sphere covering can be modeled 
by LMIs. 

In Section [8] we combine these tools and attain algorithms which theoretically solve the two problems 
under consideration. 

Due to a combinatorial explosion of the number of different Delone triangulations, our implementations 
of the algorithms only give complete solutions for d < 5. Moreover, the convex optimization algorithms 
we used are interior point methods and so yield only approximations. Therefore, in Section [9j we collect 
some mathematical tools which allow us to determine exact results from these approximations. In particular 
we can test computationally whether or not a given positive definite quadratic form gives a locally optimal 
solution. In the case of the lattice packing-covering problem, we can test if such a solution is isolated. 

In order to run a heuristic search for good lattices, it is necessary to have local lower bounds that we 
can compute fast. We describe one class of such bounds depending on the methods of inertia in Section [lOl 
which we used to find new lattices in dimension 6. 

Both problems have been previously solved only for dimensions d < 5. Our implementations not only 
verify all of these results, but also attain additional information on locally optimal solutions for d = 5. 
Moreover, we find new best known lattices for both problems in dimension 6, 7, and 8. In particular, we 
answer an open question by Ryshkov. In Section^Jwe report on our results. There we distinguish between 
conjectures for which we only have numerical evidence and theorems for which we have rigorous proofs. 
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2 Basic Concepts and Notations 



In this section we fix the notation we use throughout this paper. We refer the reader to [BR79], [GL87| and 
[CS88 ] for further information about the introduced concepts. The reader familiar with sphere packings and 
coverings, as well as lattices and their relation to positive definite quadratic forms, may skip this section. 

Let K d be a (i-dimensional Euclidean space equipped with inner product (•,•), norm || • || and unit ball 
B d = {x G E d : \\x\\ < 1}. A lattice L is a discrete subgroup in E d . From now on we assume that all 
lattices L have full rank d; that is, there exists a regular matrix A G GL^(IR) with L = A7L d . The columns 
of the matrix A are called a basis of L. All bases of L are of the form AU with U G GLd(Z). Thus, the 
determinant det(L) = | det(A)| > of the lattice L is well defined. We say that two lattices L and L' 
belong to the same isometry class if, for every basis A of L, there is a basis A' of L' and an orthogonal 
transformation O G O d (R) such that A' = OA. 

The Minkowski sum L + aB d = {v + ax : v G L, x G with a G M>o, is a lattice packing if 
the translates of aB d have mutually disjoint interiors and a lattice covering if E d = L + aB d . The pocting 
radius A(L) of a lattice L is given by 

A(L) = max{A : L + XB d is a lattice packing}, 

and the covering radius fi(L) by 

n(L) = min{/x : L + /xl^ is a lattice covering}. 

The above values are attained: The packing radius is equal to half the length of a shortest non-zero vector 
of L and the covering radius is equal to the maximum distance of points in ¥, d to a closest lattice vector. The 
packing radius is the inradius of the Dirichlet-Voronoi poly tope of L 

DV(L) = {x G E d : \\x\\ < \\x - v\\,v G L}, 

and the covering radius is its circumradius. Both functionals are homogeneous, that is, for a G M we have 

p,(aL) = \a\iJ,(L) and \(aL) = |a|A(L). 

Thus, the covering density 

is invariant with respect to scaling of L. The same is true for the packing-covering constant 

ML) 



X(LY 



All these functionals are invariants of the isometry classes. 
In this paper we study the following two problems: 



Problem 2.1. (Lattice covering problem) 

For a given d G N, determine 6^ = min^ 0(L), where L C ~E d runs over all (i-dimensional lattices. 
Problem 2.2. (Lattice packing-covering problem) 

For a given d G N, determine 7^ = min^ 'y(L), where L C K d runs over all (i-dimensional lattices. 
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We describe history and results of both problems in Section |3] and Section |4] There we assume that 
the reader is familiar with certain important lattices which are described in [CS88|, Ch. 4. Historically and 
for computational reasons, the problems were studied by using the intimate relation between lattices and 
positive definite quadratic forms (PQFs). 

Let us describe this relation: To a d-dimensional lattice L = AZ d with basis A we associate a d- 
dimensional PQF 

Q[x] = x t A t Ax = x t Gx, 

where the Gram matrix G = A 1 A is symmetric and positive definite. We abuse notation and identify 
quadratic forms with symmetric matrices by saying Q = G and Q[x] = x l Qx. The set of quadratic forms 
is a d ^- d ^ dimensional real vector space S d , in which the set of PQFs forms an open, convex cone S d . Its 
closure is the convex cone of all positive semidefinite quadratic forms S d , which is pointed at 0. 

Note that Q depends on the chosen basis A of L. Two arbitrary bases A and B of L are transformed 
into each other by a unimodular transformation, that is, there exists an U G GL^(Z) such that A = BU. 
Thus, GLrf(Z) acts on S> by Q i— > ^QU. Two PQFs lying in the same orbit under this action are called 
arithmetically equivalent. This definition naturally extends to positive semidefinite quadratic forms. 

Thus every lattice uniquely determines an arithmetical equivalence class of PQFs. On the other hand, 
every PQF Q admits a Cholesky decomposition Q = A 1 A, where the upper triangular matrix A is uniquely 
determined up to an orthogonal transformation O G O^(R). Altogether, we have a bijection between isom- 
etry classes of lattices Od(K)\GLd(M)/GI_d(Z) and arithmetical equivalence classes of PQFs 5> /GLd(Z). 

As a consequence, the lattice covering and the lattice packing-covering problem translate into problems 
for PQFs: The determinant (or discriminant) of a PQF Q is defined by det(Q). The homogeneous minimum 
\(Q) and the inhomogeneous minimum n(Q) are given by 

A(Q) = min Q[v], and n(Q) = max min Q[x — v]. 



A corresponding lattice L satisfies det(L) = det(Q), = y/^Q), A(L) = sj \{Q)/2. Therefore 
our goal is to minimize 



among all PQFs Q G S d . 

Since 0(Q) and j(Q) are invariant with respect to the action of GL^(Z) on 5> > we on ^y neec ^ to 
consider one PQF in each arithmetical equivalence class. Finding a fundamental domain in S d is one 
of the most basic and classical problems in the geometry of numbers. Such a reduction theory for PQFs, 
especially suitable for Problem IXTl and Problem IT2l is due to Voronoi. We describe it in detail in Section|5] 



3 The Lattice Covering Problem 

Kershner, in 1939, was the first to consider the lattice convering problem. In IKer3 9 1 he showed that the 
hexagonal lattice (see Figure 1) gives the most economical sphere covering in the plane even without the 
restriction of being a lattice covering. 
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Figure 1 . The sphere covering given by the hexagonal lattice. 



Since then the lattice covering problem has been solved up to dimension 5 (see Table 1). In all these 
cases the lattice A^, whose covering density equals 

e <W(l£l)^ +1) ^ 

provides the optimal lattice covering. Gameckii [Gam62|, [Gam63|, and Bleicher [Ble62| were the first to 
compute the covering density of for general d. They also showed that it is locally optimal with respect 
to covering density in every dimension. 

The optimality of the body centered cubic lattice A3 whose Dirichlet-Voronoi polytope is a regular 
truncated octahedron (the Dirichlet-Voronoi polytope of A J is a regular permutahedron) was first proven by 
Bambah [Bam54b|. Later, Barnes substantially simplified Bambah's proof in [Bar56| and strengthened the 
result by showing that in dimensions 2 and 3 the lattice k* d is the unique locally optimal lattice covering. He 
used Voronoi's reduction theory and anticipated that this is the right setup for solving the lattice covering 
problem in dimensions larger than three. Our algorithm in Section [Ol confirms his anticipation. A third 
proof, which is mainly elementary and unlike the previous two does not use any reduction theory of PQFs, 
was given by Few [F ew561 . At the moment no attempt is known to the authors to show that the optimal 
three-dimensional lattice covering also gives an optimal sphere covering without lattice restriction. 

In [Bam54a| Bambah conjectured that the lattice A4 gives the least dense four-dimensional lattice cov- 
ering. In [DR63| Delone and Ryshkov proved Bambah's conjecture. In [Bar65|, [Bar66| Baranovskii gave 
an alternative proof of this fact. He determined all locally optimal lattice coverings in dimension 4. Dickson 
[Dic67 1 gave another alternative proof of this fact. 

In a series of papers [Rys73|, [BR73|, [BR75|, [RB75| Ryshkov and Baranovskii solved the lattice 
covering problem in dimension 5. They prepared a 140-page long monograph [RB78 1 based on their inves- 
tigations. 



d 


lattice covering 


covering density 0^ 


1 


Z 1 


1 


2 


A* 


1.209199 


3 


A* 


1.463505 


4 


A* 


1.765529 


5 


A* 


2.124286 



Table 1 . Optimal lattice coverings. 
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In [ Rys67 1 Ryshkov raised the question of finding the lowest dimension d for which there is a better 
lattice covering than the one given by k* d . In the same paper he showed that fk* d is not the most efficient 
lattice covering for all even d > 114 and for all odd d > 201. One of our main results in this paper is the 
answer to Ryshkov 's question (see Theorem II 1.3b : Dimension d = 6 is the lowest dimension for which 
there is a better lattice covering than the one given by A^. The proof is based on a computer search. In 
Section HTTl we will give an algorithm which finds all locally optimal lattice coverings in a given dimension. 
Using this algorithm we are able to verify all of the known results about optimal lattice coverings up to 
dimension 5. Unfortunately, due to a combinatorial explosion, the algorithm can not be applied practically 
in dimension 6 or greater. Nevertheless, we are able to find good lattice coverings in dimension 6, 7, and 8 
by applying several heuristics. We will give more details in Section fTTI 



d 


lattice 


covering density 


6 


icl 
L 6 


2.464801 


7 


1 c 
L 7 


2.900024 


8 


LI 


3.142202 


9 


A| 


4.340185 


10 


A* 


5.251713 


11 


A 4 


5.598338 


12 


A* 
M 12 


7.510113 


13 


A 7 
M 13 


7.864060 


14 


A 5 


9.006610 


15 


A 8 


11.601626 


16 


A* 
M 16 


15.310927 


17 


A* 
M 17 


18.287811 


18 


A* 
M 18 


21.840949 


19 


A* 


26.081820 


20 


A* 


31.143448 


21 


A* 


37.184568 


22 


A* 
A 22 


< 27.8839 


23 


A* 
A 23 


< 15.3218 


24 


A24 


7.903536 



Table 2. Least dense known (lattice) coverings up to dimension 24. 

What else is known? In Table 2 we list all the least dense known lattice coverings in dimensions 6 
to 24. At the same time this list gives the least dense known sphere coverings: There is no covering of equal 
spheres known which is less dense than the best known lattice covering. This table is an update of Table 2.1 
in [CS88] and we provide an up-to-date table on our web page [SV05]. We conclude this section by briefly 
describing the origins of these updates. 

The Leech lattice A24 yields the best known lattice covering in dimension 24. The covering density 
of the Leech lattice was computed by Conway, Parker and Sloane ([CS88| Ch. 23). It is not too brave to 
conjecture that the Leech lattice gives the optimal 24-dimensional sphere covering. In [SV04] we took a 
first step towards proving this conjecture by showing that the Leech lattice gives a locally optimal lattice 
covering. Using the Leech lattice Bambah and Sloane constructed in [BS82 1 a series of lattices in dimensions 
d > 24 which give a thinner lattice covering than AJ. It seems that as a "corollary" of the existence of the 
Leech lattice the duals of the laminated lattices A22 and A23 give good lattice coverings. Their covering 
densities were estimated by Smith [Smi88|, but we do not know the exact values of 0(A| 2 ) and ©(A^). 
For the definitions of the root lattices A n , D n , E n , and the laminated lattices A n , we refer to [CS88 1, Ch. 4 
and Ch. 6. 
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In ICox51l Coxeter gave a list of locally optimal lattice packings (extreme lattices) which are related to 
Lie groups. One of them is the infinite series of locally densest packing lattices where d > 2 and r > 1 
divides d + 1. The lattice is the unique sublattice of containing to index r. In [Bar94] Baranovskii 
determined the covering density of the lattice covering given by Ag, which is slightly better than the one 
given by Ag. Recently, Anzin extended Baranovskii 's work. In [ Anz02] he computed the covering densities 
of A^ and A^ 3 , and in a private communication he reported on computing the covering densities of Af 4 and 
Af 5 . They all give less dense lattice coverings than those provided by the corresponding A^. We do not 
know whether these lattice coverings are locally optimal. 

To answer Ryshkov's question exhaustively it will be necessary to further investigate lattice coverings in 
the dimensions d = 10, 12, 16, 17, ... , 21. We hope that the methods we present in this paper will be useful 
for this project. 

4 The Lattice Packing- Covering Problem 

The lattice packing-covering problem has been studied in different contexts and there are several different 
names and interpretations of the lattice packing-covering constant 7^. Lagarias and Pleasants [LP02] re- 
ferred to it as the "Delone packing-covering constant". Ryshkov |Rys74| studied the equivalent problem of 
minimizing the density of (r, i?)-systems. An (r, i?)-system is a discrete point set X C M. d where (1) the 
distance between any two points of X is at least r and (2) the distance from any point in M. d to a point in X 
is at most R. If X is a lattice, then r = X(X) /2 and R = fi(X). 



Figure 2. A close sphere packing given by the hexagonal lattice. 

Geometrically we may think of solving the lattice packing-covering problem as of maximizing the min- 
imum distance between lattice points in a lattice covering with unit spheres. Alternatively, we may think 
of it as minimizing the radius of a largest sphere that could additionally be packed into a lattice packing 
of unit spheres (see Figure 2). This minimal gap-radius is equal to 7^ — 1. Therefore the problem raised 
by L. Fejes Toth [T6t76| of finding "close packings" attaining this gap-radius is another formulation of the 
packing-covering problem. 

The last interpretation shows that 7^ > 2 would imply that in any d-dimensional lattice packing with 
spheres of unit radius there is still space for spheres of radius 1. In particular, this would prove that densest 
sphere packings in dimension d are non-lattice packings. This phenomenon is likely to be true for large 
dimensions, but it has not been verified for any d so far. 

Problem 4.1. Does there exist a d such that 7^ > 2? 
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Note that this problem is particularly challenging in view of the asymptotic bound 7^ < 2 + o(l) due to 
Butler IBut72l . 

As for the lattice covering problem, the lattice packing-covering problem has been solved up to dimen- 



sion 5 (see Table 3). Ryshkov [Rys74] solved the general 2-dimensional case. The 3-dimensional case was 
settled by Boroczky [Bor86], even without the restriction to lattices. The 4- and 5-dimensional cases were 
solved by Horvath [Hor82|, [Hor86|. Note that the lattices H04 and H05 (see Section ITU discovered by 
Horvath are neither best covering, nor best packing lattices. As in the case of lattice coverings, the results 
were attained by using Voronoi's reduction theory. 



d lattice lattice packing-covering constant 7^ 



d 


lattice 


1 


Z 1 


2 


A* 
A 2 


3 


A* 


4 


H04 


5 


H05 



1 

4? « 1.154700 

v3 



5/3 « 1.290994 
V gV^yg - 1) ss 1.362500 
3/2 + v/13/6 ~ 1.449456 



Table 3. Optimal packing-covering lattices. 

Our computations, described in Section^TJ verify all of the known results in dimension < 5. As for the 
covering problem, none of the values 7^ has been determined in a dimension d > 6 so far. In Section [TT 
we report on a new best known packing-covering lattice for d = 6. We thereby show in particular that 76 < 
1.412, revealing the phenomenon 76 < 75, recently suspected by Lagarias and Pleasants [LP02], Sec. 7. 
Note that this was already observed by Zong [Zon02] (Remark 3), who showed that 7(Eg) = \/2 < 75. 

We were not able yet to find any new best known lattices in dimensions d > 7. The lattice E£ gives 
the best known lattice in dimension 7. Nevertheless, because of their symmetry and the known bounds on 
7d, Zong [Zon02| (Conjecture 3.1) made the following conjectures: Eg and Leech lattice A24 are optimal 
in their dimensions. In IS V04II we showed that the Leech lattice gives a locally optimal lattice packing- 
covering constant. It is an open question whether or not the root lattice Eg gives a locally optimal lattice 
packing-covering constant as well. The corresponding values of the lattice packing-covering constant are 
shown in Table 4. In dimension 7 and between dimensions 8 and 24 we do not yet know enough to state any 
serious conjectures. The exact value for the smallest known lattice packing-covering constant in dimension 6 



is 2V2^ 



56 (see SectionlTTli. 



d 


lattice 


lattice packing-covering constant 7^ 


6 




2v / 2\/798 - 56 ss 1.411081 


7 


E7 


y/Y]?> ss 1.527525 


8 


E 8 


1.414213 


24 


A24 


y/2 « 1.414213 



Table 4. Lattice packing-covering records. 



5 Voronoi's Reduction Theory 

The general task of a reduction theory for PQFs is to give a fundamental domain for cS^q/GL^Z). This is a 
subset which behaves like 5> /GL^(Z) up to boundary identifications. There are many different reduction 
theories, connected with names like Lagrange, GauB, Hermite, Korkine, Zolotareff, Minkowski, Voronoi, 
and others (see [S085]). In this section we describe the reduction theory developed by Voronoi in [Vor08|. 
It is based on Delone triangulations. 
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5.1 Secondary Cones of Delone Triangulations 

Let Q £ S d be a PQF. A polytope L = convj^i, . . . , v n }, with V\,...,v n in Z d , is called a Delone 
poly tope of Q if there exists acei d and a real number r£l with Q[vi — c] = r 2 for alH = 1, . . . , n, 
and for all other lattice points v £ . . . , i? n } we have strict inequality Q[v — c] > r 2 . The set of all 

Delone polytopes 

Del(Q) = {L : L is a Delone polytope of Q} 

is called the Delone subdivision (or L-partition) of Q. A Delone triangulation is a Delone subdivision that 
consists of simplices only. For historical reasons we use the letter L to denote Delone polytopes (see [ Vor08 1 
and IIDeDHn . 

The Delone subdivision of a PQF is a periodic polytopal subdivision of M. d . We say that two Delone 
polytopes L and L' are equivalent if there is a v £ Z rf with L' = L + v. Given a Delone subdivision V of 
IR rf , the set of PQFs with Delone subdivision V forms the secondary cone 

A(V)={Q £S d :Del(Q)=V}. 

In the literature it is often referred to as the L-type domain of Q £ A (2?). We prefer the term "secondary 
cone" because of the close connection of Voronoi's reduction theory to the theory of secondary polytopes 
which we will point out in Section 1531 

Let Q be a PQF whose Delone subdivision is a triangulation of M. d . In the following we will describe 
the secondary cone of Del(Q). For this, let L = conv{t>i, . . . ,Vd + i} and L' = conv{t>2, • • • , 
be two d-dimensional Delone simplices of Q sharing the common facet F = conv{?;2, • • • Let 
ai,..., a>d+2 be real numbers with ct\ = 1, Yli=i «i = and J2i=i a i v i = (hence ctd+2 > 0). The 
regulator of the pair of adjacent simplices (L, L') is the linear form Q(l,l') (Q') = Yli=i a iQ'[ v i\> Q' £ S d . 
In particular, note that the regulator solely depends on the points v\, . . . ,Vd+2, that Q(l,u){Q) > 0, and 
that Q(l,L') = Q(L+v,L'+v) f° r ai l v ^ <^ d - One can describe A(Del(Q)) by linear inequalities coming from 
the (finitely many) regulators of Del(Q): 

Proposition 5.1. ( IIVor0 8 1, §77) Let Q be a PQF whose Delone subdivision is a triangulation. The secondary 
cone of the Delone triangulation Del(Q) is the full-dimensional open polyhedral cone 

A(Del(Q)) = {Q' £ S d : Q^ L >){Q') > 0, for all pairs (L, V) of adjacent simplices of Del(Q)}. 

Note on the one side that A(Del(Q)) is contained in 5> by definition and on the other side that 
the linear inequalities given by the regulators imply that a quadratic form which satisfies them is positive 
definite. 

Example 5.2. As a first example and because of its importance for the lattice problems introduced in Sec- 
tion |2j we describe the Delone subdivision of Voronoi's principal form of the first type Q[x] = dJ2 x f ~ 
Y^XiXj, which is associated to the lattice A^, in greater detail. The Delone subdivision of Q is a tri- 
angulation and can be described as follows: Let e±, . . . , be the standard basis vectors of Z d , and set 
e^+i = — ei — • • • — e<f. For a permutation it £ S^+i we define the d-dimensional simplex L v by 

L n = cony{e n ^, + e n(2) , . . . , e A1) -\ h e^^}. 

The set of simplices {L w + v : v £ Z d , tt £ S^+i} defines a triangulation of M. d which we from now on 
denote by T>\. The full-dimensional cells containing the origin are L w , tt £ S^+i- Two simplices L n and L^i 
have a facet in common if and only if n and n' differ by a single transposition of two adjacent positions. The 
automorphism group of T>\ is isomorphic to the permutation group S^+i. The star of the origin is illustrated 
in Figure 3. This consists of all Delone polytopes containing the origin. 
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The secondary cone of T>\ is (see IVor08l . §102-104) 

A(Di) = {Q £ S d : qij < for i ^ J and > for j = 1, . . . , d }. 

Its topological closure A(2?i) is called Voronoi's principal domain of the first type. □ 

It was shown by Voronoi ([Vor08], §97) that the topological closures of the secondary cones gives a 
facet-to-facet tessellation of S> . By a theorem of Gruber and Ryshkov [GR89| we even have a face-to-face 
tessellation because "facet-to-facet implies face-to-face". This means that a face sharing relative interior 
points with another face of its dimension coincides with this face, whenever this property holds for the 
facets (faces of co-dimension 1). 

5.2 Bistellar Operations 

Now, given a secondary cone, how do we find its neighbors, that is, those secondary cones sharing a facet 
with the given one? 

An answer can be given by taking a closer look at the so-called repartitioning polytopes (introduced by 
Ryshkov and Baranovskii in [RB78 1), which are "hidden" in the definition of the regulators. Repartitioning 
polytopes are d-dimensional Delone polytopes having a representation as the convex hull of two Delone 
simplices sharing a common facet. Thus, repartitioning polytopes have d + 2 vertices. 

Generally, d-dimensional polytopes with d + 2 vertices have the special property that there are ex- 
actly two different ways to triangulate them: Let V be a set of d + 2 points which affinely spans IR d . Let 
Yl,v&v OLvV = 0> J2 v ev a v = 0' be an affine relation between these points. There exist exactly two tri- 
angulations of convV: T + (V,a) with (i-simplices conv(V\{«}), a v > 0, and 71 (V, a) with d-simplices 
conv(T/\{t?}), a v < (see for example [GKZ94|, Ch. 7, Prop. 1.2). 

Let V be a Delone triangulation of M. d and let F be a (d— 1) -dimensional cell of V. Then, F is contained 
in two simplices L and V of V. By V we denote the set of vertices of L and L',V = vert L U vert L'. By 
oc we denote an affine relation between the points in V. The (d — 1) -dimensional cell F is called aflippable 
facet of the triangulation T> if one of the triangulations T + (V, a) or 71 (V, a) is a subcomplex of V. If F 
is a flippable facet of V and we replace the subcomplex T + (V,ct) by 71 (V, a) [respectively 71 (V, a) by 
1+(V, a)], then we get a new triangulation. This replacement is called bistellar operation or flip. 

Notice that non-flippable facets do exist and that performing a bistellar operation in a Delone triangula- 
tion does not necessarily produce a Delone triangulation. Both phenomena occur starting from dimension 4. 
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Figure 4. The graph of two-dimensional Delone triangulations. 



Nevertheless, the facets of A(P) correspond exactly to those bistellar operations of V which yield new 
Delone triangulations. A (d— 1) -dimensional cell LnL' G T> is a flippable facet whenever the corresponding 
regulator Q(l,l') gives a facet-defining hyperplane of A(D) (see IVor08l . §87-88). This is clear since the 
repartitioning polytope conv(L U L) is a Delone polytope of the PQFs lying in the relative interior of the 
facet given by Q^vy 

Let F be a facet of A(P). We describe how the Delone triangulation V changes if we vary a PQF 
continuously: We start from the interior of A(£>), then we move towards a relative interior point of F and 
finally we go infinitesimally further, leaving A(D). In every repartitioning polytope V = conv(L U L r ) 
where L,L' is a pair of adjacent simplices whose regulator defines F, i.e. the linear span of F satisfies 
linF = {Q G S d : Q(l,L'){Q) = 0}, we perform a bistellar operation. This gives a new triangulation V, 
which is again a Delone triangulation. The two secondary cones A(D) and A(D') have the complete 
facet F in common. We say that V and V are bistellar neighbors. In [Vor08] §91-96, Voronoi computes 
the secondary cone of V explicitly and shows that A(D') has dimension d i d ^' 1 > , 

5.3 Main Theorem of Voronoi's Reduction Theory 

By constructing bistellar neighbors we could produce infinitely many Delone triangulations starting from 
the Delone triangulation T>\ of Voronoi's principal form of the first type (a part of the infinite flip graph 
of two-dimensional Delone triangulations is given in Figure 4). Many of these will not be essentially new, 
because the group GL^(Z) is acting on the set of Delone subdivisions by (A,V) i— ► AD and it is acting on 
the set of secondary cones by (A, A) h-> A l A.A for A G GL^(Z). We are only interested in the orbits 
of these group actions and there are only finitely many, as shown by Voronoi [Vor08| §98 (see also Deza, 
Grishukhin and Laurent in [DL97|, Chapter 13.3). Altogether we get: 

Theorem 5.3. {Main Theorem of Voronoi's Reduction Theory) 

The secondary cone of a Delone triangulation is a full-dimensional, open polyhedral cone in S d . The 
topological closures of secondary cones of Delone triangulations give a face-to-face tiling of S d . Two 
secondary cones share a facet if and only if they are bistellar neighbors. The group GL^(Z) acts on the 
tiling, and under this group action there are only finitely many inequivalent secondary cones. 

The main theorem translates into Algorithm 1 which enumerates all inequivalent Delone triangulations 
in a given dimension. We developed the program sec (secondary cone cruiser) which is an implemen- 
tation of Algorithm 1. The interested reader can download sec from our web page ISV05I . Using our 
implementation we succeeded in reproducing the classification of all inequivalent Delone triangulations up 
to dimension 5. Table 5 shows the numbers. 
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Input: Dimension d. 

Output: Set 1Z of all inequivalent d-dimensional Delone triangulations. 
T <— {T>{\, where T>\ is the Delone triangulation described in Example l5 .2! 

while there is a V G T do 

T «- 7e «- ft U {£>}. 

compute the regulators of V. 
compute the facets Fx, . . . ,F n of A(X>). 
for i = 1, . . . , n do 

compute the bistellar neighbor T>i of V which is defined by F{. 

if T>i is not equivalent to a Delone triangulation in 1Z U {Pi, . . . , Pj-i} then 

r<-Tu{Pj}. 

end if 
end for 

end while 

Algorithm 1 . Enumeration of all inequivalent Delone triangulations. 



We are not the first to compute this classification. Voronoi performed the classification of all inequiv- 
alent Delone triangulation up to dimension 4 in his memoir [Vor08]. In [BR73] and [RB78] Ryshkov and 
Baranovskii reported on 221 inequivalent Delone triangulation in dimension 5. However, they missed one 
type which was found by Engel [Eng98|. Grishukhin and Engel [EG02] undertook the non-trivial task of 
identifying the Delone triangulation missing in the list of Ryshkov and Baranovskii. There they also re- 
port on several errors in both lists. Our computations confirm that the number of Delone triangulations in 
dimension 5 is 222. Beginning with dimension 6 the number of inequivalent Delone triangulations starts 
to explode. We found more than 250, 000 inequivalent triangulations before our implementation stopped 
because of memory reasons. We do not know how many inequivalent triangulations we have to expect in 
dimension 6 but we are quite certain that we only saw a small fraction of them. No non-trivial bounds on 
the number of inequivalent Delone triangulations for general dimension d are known. 



Dimension 


1 


2 


3 


4 


5 


# Delone triangulations 


1 


1 


1 


3 


222 



Table 5. Numbers of inequivalent Delone triangulations. 



5.4 Degeneracy 

Until now we have only dealt with Delone triangulations of positive definite quadratic forms. Let us look at 
possible degenerations — Delone subdivisions of positive semidefinite quadratic forms — and find out how 
they fit into the theory developed so far. 

Let Q be a positive semidefinite quadratic form which is arithmetically equivalent to ( o o' ) where Q' 
is positive definite. Then, we can use the definition of Delone subdivision almost literally; we only have to 
replace "polytope" by "polyhedron" (a polyhedron is the intersection of finitely many half spaces; apolytope 
is a bounded polyhedron). 

Delone subdivisions are limiting cases of triangulations. Their secondary cones occur on the boundaries 
of full-dimensional secondary cones of Delone triangulations. Let V and V be two Delone subdivisions. 
We say V is a refinement of V if every Delone polytope of V is a subset of some Delone polytope of V. The 
following proposition, which seems to be folklore, shows that the relation between refinements, secondary 
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cones and sums of positive semidefinite quadratic forms is very natural. One can find a proof for example 
in IVaKBl fProo. 2.6.1). 



Proposition 5.4. Let V be a Delone triangulation. 

1. A positive semidefinite quadratic form Q lies in A(D) if and only if I? is a refinement of Del(Q). 

2. If two positive semidefinite quadratic forms Q and Q' both lie in A(P), then Del(Q + Q') is a 
common refinement of Del(Q) and Del(Q'). 

Therefore, the classification of all inequivalent Delone subdivisions is equivalent to the classification of 
all inequivalent secondary cones. This has been done up to dimension 5. The 1- and 2-dimensional cases are 
trivial; the 3-dimensional case goes back to Federov, who classified all polytopes which tile 3-dimensional 
space by translates in 1885. Delone [Del29| (later corrected by Stogrin [Sto73|) found 51 of the 52 Delone 
subdivisions in dimension 4. Recently, Engel [EngOO] reported that there are 179,372 inequivalent five- 
dimensional Delone subdivisions. It is possible to verify his result by enumerating the faces of the 222 
secondary cones of Delone triangulations in dimension 5 up to equivalence. Since we do not need Engel's 
result for our application we did not verify it. Again, the number in dimension 6 is not known and will be 
much larger. Table 6 summarizes this discussion. 



Dimension 


1 


2 


3 


4 


5 


# Delone subdivisions 


1 


2 


5 


52 


179,372 



Table 6. Numbers of inequivalent Delone subdivisions. 



5.5 Generalized Secondary Polytopes 

Triangulations of discrete point sets have attracted many researchers in recent years. They have many 
applications, for example in computational geometry, optimization, algebraic geometry, topology, etc. One 
main tool to understand the structural behavior of triangulations of finite point sets is the theory of secondary 
polytopes invented by Gel'fand, Kapranov and Zelevinsky ([GKZ94], Ch. 7). 

Let A = {oi, . . . , a n } C R d be a finite set of points. Let w : A — > R be a map that assigns to 
every point in A a weight. The set of weight maps forms a vector space over R which we denote by R- 4 . 
A lifting map I : A -> R d x R, Z(a») = (ai,w(a,i)) is defined by w which lifts each point etj £ A 
by its weight w(a,i). A subdivision of the convex polytope conv .4. is induced by I: We take the convex 
hull of the lifted points conv 1(A) and project its lower faces as seen from (0, — oo) back down onto R d . A 
subdivision that can be obtained in this manner is called a regular subdivision. Delone subdivisions (or more 
precisely Delone subdivisions of finitely many points) are regular subdivisions since the underlying positive 
semidefinite quadratic form is the weight function. This view on Delone subdivisions was introduced by 
Brown [Bro79] and by Edelsbrunner and Seidel [ES86|. 

Let T be a regular triangulation of conv A. We may ask what are the weight functions which define T. 
What is the secondary cone of T in the parameter space R-^? As in Voronoi's reduction theory it turns out 
that the secondary cone of T is a full-dimensional open polyhedral cone. The topological closures of the 
secondary cones of all regular triangulations tile the space R- 4 face-to-face. The tiling is called secondary 
fan of A. If two secondary cones have a facet in common, then the corresponding regular triangulations 
differ by a bistellar operation in exactly one "repartitioning polytope" (in this context it is a polytope with 
d + 2 vertices without the condition of being a Delone polytope) that is defined by the facet. The faces in 
the secondary fan A are in a one-to-one correspondence to regular subdivisions in essentially the same way 
we discussed in Section l5~4l for Delone subdivisions. 
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So far we have seen that the theory of regular subdivisions of finite point sets and the theory of Delone 
subdivisions of the lattice 7L d can be developed analogously, but there are also differences: The parame- 
ter space for regular subdivisions is the vector space R" 4 , while for Delone subdivisions it is the pointed 
cone S d . Groups play an important role for Delone subdivisions. The group Z d is acting on Delone 
subdivisions by translations. On the set of secondary cones the group GL^(Z) is acting. 

If we order all regular subdivisions of conv A by refinement we get a poset. This poset has a very nice 
combinatorially structure as proved by Gel 'f and, Kapranov and Zelevinsky : There exists a poly tope — the 
secondary poly tope S(.A) of A — whose normal fan equals the secondary fan of A. So the refinement poset 
is anti-isomorphic to the face lattice of the secondary polytope. Regular triangulations are in one-to-one 
correspondence to the vertices, two regular triangulations differ by a bistellar operation if and only if their 
vertices are connected by an edge, etc. 

Recently, a similar combinatorial structure lurking behind the refinement poset of Z d -periodic subdivi- 
sions (this poset contains the poset of Delone subdivisions) has been described by Alexeev [ Ale02], Sec. 5.11 
and 5.12 . He gives an unbalanced and a balanced version of these generalized secondary polytopes, where 
the latter one is invariant with respect to the group action of GI_d(Z). 



6 Convex Optimization with LMI Constraints 

In this section we introduce determinant maximization problems, which are convex programming problems 
with linear matrix inequality constraints. In a sense they are equivalent to the better known semidefinite pro- 
gramming problems. For both classes efficient algorithms and implementations are available. In Section [8] 
we will see how Problem |2~T1 and Problem [2.2l can naturally be formulated as a finite number of determinant 
maximization problems. 

Following Vandenberghe, Boyd, and Wu [VBW98| we say that a determinant maximization problem is 
an optimization problem of the form 



minimize c l x — log det G(x) 
subject to G(x) is a positive definite matrix, 

F(x) is a positive semidefinite matrix. 



(1) 



The optimization vector is x G 
R nxn are affme maps: 



where G -, G 



and Fs G 



l d , the objective vector is c G R d , and G : M d -> R mxm and F : R d -» 

G(x) = G + x x Gi + • • • + x d G d , 
F(x) = F + xiF 1 + --- + x d F d , 

xn , for i = 0, . . . , d, are symmetric matrices. In the sequel we will write 



G(x) y and F(x) y for the linear matrix inequalities defining the constraints of the determinant 
maximization problem. As a special case, our formulation reduces to a semidefinite programming problem 
whenever G(x) is the identity matrix for all x G R rf . 

Currently there exist two different types of algorithms — ellipsoid and interior-point methods — which 
efficiently solve semidefinite programming problems. They can approximate the solution of a semidefinite 
programming problem within any specified accuracy and run in polynomial time if the instances are "well- 
behaved". For more information on the exciting topic of semidefinite programming the interested reader is 
referred to the vast amount of literature which to a great extend is available on the World Wide Web. A good 
starting point is the web page 1 of Christoph Helmberg. 



^http : / / www-user . tu-chemnit z . de/ "helmberg/ semidef ■ html| 
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Nesterov and Nemirovskii [ NN94 ] developed a framework for the design of efficient interior-point algo- 
rithms for general and specific classes of convex programming problems. There (§6.4.3), they also showed 
that the determinant maximization problem can be transformed into a semidefinite programming problem 
by a transformation which can be computed in polynomial time. Nevertheless it is faster to solve the 
determinant maximization problem directly. Vandenberghe, Boyd, Wu [VBW98| and independently Toh 



[Toh99| gave interior-point algorithms for the determinant maximization problem. Both algorithms fit into 
the general framework of Nesterov and Nemirovskii. For our implementation we use the software package 
MAXDET of Wu, Vandenberghe, and Boyd as a subroutine. 

One nice feature of determinant maximization problems is that there is a duality theory similar to the 
one of linear programming, which allows one to compute certificates for a range in which the optimum is 
attained, the so called duality gap. The dual problem of problem Q is (see [ VBW98 1, Sec. 3) 



maximize log det W — trace(GoW) — trace(FoZ) + m 
subject to ti&ce(GiW) + trace(FjZ) = Cj for i = 1, . . . , d, 

w y o, z y o. 



(2) 



Here, W £ ]R mxm and Z G M nxn are symmetric matrices of the same size as the G-i and Fi respectively. 

The knowledge of a vector x with G(x) y and F(x) y 0, and of a pair (W, Z), with trace(GjVF) + 
trace(FjZ) = Cj for i = 1, . . . , d, and W y 0, and Z y 0, gives upper and lower bounds for the optimal 
value p* of by 

log det W - trace(G TV) - trace(F Z) + m < p* < c l x - logdetG(x). 

If x, W and Z have rational entries only, we can compute lower and upper bounds of p* that are mathemati- 
cal rigorous. For the proofs in Section[^we have therefore developed a program rmd (rigorous MAXDET) 
which calls MAXDET to find floating point values of a vector x and a pair (W, Z) satisfying the constraints. 
Then it approximates x and (W, Z) by rational numbers and checks whether these rational approximations 
satisfy the constraints. In such a case they guarantee a certified duality gap. The interested reader can 
download rmd from our web page I SV05I . 



7 An LMI Constraint for the Inhomogeneous Minimum 

In this section we will give a linear matrix inequality in the parameters (qij) of a PQF Q = (qij), which is 
satisfied if and only if the inhomogeneous minimum of Q is bounded by a constant, say p,{Q) < 1. For this 
it is crucial to observe that x t Qy is a linear expression in the parameters (qij). The PQF Q = (q^) defines 
the inner product of a Euclidean space (M. d , (-,■)) by (x, y) = x t Qy. 

From the article [DDRS70|, §3, of Delone, Dolbilin, Ryshkov and Stogrin we can extract the following 
proposition which will be central in our further discussion. 

Proposition 7.1. Let L = conv{0, v\, . . . , Vd} Q M. d be a d-dimensional simplex. Then the radius of the 
circumsphere of L is at most 1 with respect to (•, •) if and only if the following linear matrix inequality (in 
the parameters gy) is satisfied: 





f 4 


Oi, 


Vi) 


(V2,V 2 ) ■ 


■ (Vd,V d )\ 




(v 1 ,v 1 ) 


{vi, 


Vi) 


{vi,v 2 ) . 


■ (vi,v d ) 


BRl(Q) = 




(t>2, 


Vi) 


(v 2 ,v 2 ) ■ 


• («2,«d) 




\{Vd,V d ) 




Vx) 


(v d ,v 2 ) ■ 


• (v d ,v d )/ 



^http: //www. Stanford. edu/~boyd/MAXDET.html| 
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In [DDRS70| Delone et al. used this proposition as the key ingredient for showing that the set of PQFs 
which determine a circumsphere of L with radius of at most R is convex and bounded (see Proposition l9.lt . 
Because of the importance of the proposition and for the convenience of the reader we give a short proof 
here. 

Proof. We make use of Cayley-Menger determinants. The Cayley-Menger determinant of d + 1 points 
vq, . . . , v d with pairwise distances dist(vi,Vj) is 



CMOo, ...,v d ) 



1 

1 distOo,^o) 2 



dist(v ,«d) 



1 dist(> d ,i>o) 2 ... dist(v d: v d ) 2 
The squared circumradius R 2 of the simplex L equals (see for example [Ber87 1 Prop. 9.7.3.7) 

I det(dist(^,^) 2 ) < ij < d 



CM(v Q ,...,v d ) 



Replacing dist(aj, y) 2 by (x, x) — 2(x, y) + (y, y), using vq = 0, and performing elementary transforma- 
tions of the determinants turns the above formula into 

(«2,«2) ••• (Vd,V d ) 

(vi,vi) (v 1 ,v 2 ) ... (vi,Vd) 



R 2 



(v d ,v d ) (v d ,vi) (v d ,v 2 ) ... (v d ,v d ) 



The inequality R < 1 is equivalent to 



4-det((i; i ,« i )) 1<i , <d + 



det ((»i»^))l<i,,-<d 



(«2,«2) 
(«1,«2) 



(3) 



(«d,«d) 



{v d ,v d ) (v d ,v 1 ) (v d ,v 2 ) ... («d,«d) 



>0, 



which is equivalent to 



4 Oi,fi) (v 2 ,v 2 ) 



(v d ,v d ) 



(v 1 ,v 1 ) (v 1 ,v 1 ) {v 1 ,v 2 ) ... (vi,v d ) 



(v d ,v d ) (v d ,vi) (v d ,v 2 ) ... (v d ,v d ) 



> 0. 



The minors in the lower right are all determinants of Gram matrices and therefore non-negative. Hence, the 
matrix BR^(Q) is positive semidefimte. □ 

Corollary 7.2. For any d-dimensional simplex L C M. d with vertex at and Q G 5> we have 



\BR L (Q)\ >o 



BRl(Q) h 0. 
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Example 7.3. Let us compute the matrix linear inequality BR^(Q) h 0, Q = {qll qll), for the two- 
dimensional simplex L = conv{ ( Q ) , (*) , (J) }. We have 



4 gn gn + 2q 2 i + 922^ 

BRl(Q) = ( 9n 9n 911 + 921 

K q u + 2g 2 i + 922 9n + 921 9n + 2921 + 922, 

'4 0\ /0 1 1\ /0 2 N 

Ugii 1 1 1 +92i 1 

,0 0/ \1 1 1/ \2 1 2/ 



+ 922 




□ 



For a Delone polytope L other than a simplex the circumradius is less than or equal to 1 if and only if 
it is less than or equal to 1 for some d-dimensional simplex V with vertices in vert L. Therefore we 
set BRl(Q) = BKl'(Q) in this case. Since a block matrix is semidefinite if and only if the blocks are 
semidefinite, we have the following proposition which allows us to express the constraint "//(Q) < 1" 
locally by a single LMI. 

€ 5> be a PQF and let Li, . . . , L n be a representative system of all 



Proposition 7.4. Let Q = (qij 

inequivalent <i-dimensional Delone polytopes in Del(Q). Then 



BR Ll (Q) 







/i(Q) < 1 



BRl 2 (Q) 








>- 0. 



BRl„(Q) / 



8 Algorithms 

In this section we present algorithms which in theory solve the lattice covering problem and the lattice 
packing-covering problem in any dimension d. 

Our algorithm for the lattice covering problem computes all locally optimal lattice coverings of a given 
dimension. These are only finitely many because we will see that for every fixed Delone triangulation V 
there exists at most one PQF lying in the topological closure of the secondary cone of V giving a locally 
optimal covering density. So, we fix a Delone triangulation and try to find the PQF which minimizes the 
density function in the topological closure of the secondary cone of the fixed Delone triangulation. We will 
formulate this restricted lattice covering problem as a determinant maximization problem. 

Our algorithm for the lattice packing-covering problem operates similarly. For every Delone triangula- 
tion we have to solve a semidefinite programming problem. 



8.1 Solving the Lattice Covering Problem 

Recall that the covering density of a PQF Q in d variables is given by @(Q) = ^ ^}®q • ^d- Scaling of Q 
by a positive real number a leaves 6 invariant. Consequently, we can restrict our attention to those PQFs Q 
with n(Q) = 1. Hence, we solve the lattice covering problem if we solve the optimization problem 



maximize 


det(Q) 


subject to 


Q is a positive definite matrix, 




KQ) = i, 
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where the optimization variables qij are the entries of the PQF Q. The major disadvantage of this optimiza- 
tion problem is that the second constraint is not expressible as a convex condition in the variables and 
that the problem has many local maxima. A locally optimal solution is also called a locally optimal lattice 
covering. 

We will circumvent this by splitting the original problem into a finite number of determinant maximiza- 
tion problems. For every Delone triangulation T> we solve the optimization problem 



maximize 


det(Q) 


subject to 


QeA(p), 




KQ) < i- 



The relaxation of no longer requiring p,(Q) = 1 in the third constraint does not give more optimal solutions 
because if Q satisfies the constraints, then so does j^q^Q- Now, we have to show that this is indeed a 
determinant maximization problem. We have seen in Proposition 15 . 1 I that the first constraint can be expressed 
with inequalities linear in q^. The constraint fi(Q) < 1 can be expressed by a linear matrix inequality as in 
Proposition 17.41 

The optimization vector x e W d ( d+1 ^ 2 is the vector of coefficients of Q. The linear matrix inequality 
G(x) y is given by G{x) = Q. We encode the two other constraints Q G A(Q) and /i(Q) < 1 by block 
matrices in the linear matrix inequality ^(a;) t: 0. For any linear inequality which is needed to describe the 
secondary cone we have a 1 x 1 block matrix. For any inequivalent (i-dimensional simplex L G V we have 
the (d + 1) X {d + 1) block matrix BR L (Q). 

8.2 Solving the Lattice Packing- Covering Problem 

Along the same lines as above we formulate the lattice packing-covering problem as a finite number of 
semidefmite programming problems. 

Recall that the packing-covering constant of a PQF Q is jiQ) = 2 • J j^j- Since 7 is homogeneous 
we can again assume fi(Q) = 1, and the lattice packing-covering problem is equivalent to the following 
optimization problem. 



maximize 


HQ) 


subject to 


Q is a positive definite matrix, 




KQ) = i- 



A locally optimal solution of the optimization problem is called a locally optimal lattice packing-covering. 

From the previous discussion we know how to deal with the constraint n(Q) = 1. Now, how do we 
maximize A(Q)? We say that v £ Z d \{0} is a shortest vector of Q if Q[v] = X(Q). A theorem of Voronoi 
(see IVor0 8 1 §55, or [CS8SJ Ch. 21, Th. 10) implies that a shortest vector v gives the edge [0,v] in the 
Delone subdivision of Q. In a fixed Delone subdivision V there are only finitely many (at most 2(2 d — 1), 
see IVor08l §55) edges of the form [0, v]. We can maximize X(Q) for Q G A(P) as follows: We introduce 
a new variable C which we want to maximize subject to the constraints Q[v] > C where v runs through all 
edges of the form [0, v] in V. This assures X(Q) = C when C attains a maximum. All these expressions are 
linear in the coefficients qij of Q. Hence we have to solve the following semidefmite programming problem 
for every inequivalent Delone triangulation V in order to solve the lattice packing-covering problem: 
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maximize 


C 


subject to 


Q € A(V), 




KQ) < i. 




Q[v] > C, where [0, v] is an edge in V. 



9 Local Optima 

Since we are dealing with convex optimization problems, we can extract some structural information about 
uniqueness and invariance properties of locally optimal solutions. On the one hand these help us to identify 
exact coordinates of optimal solutions. On the other hand they allow us to decide whether we have found an 
isolated locally optimal solution or not. 



9.1 Properties of Local Optima 

Most of the desired properties follow from the special structure of sets with PQFs Q attaining certain fixed 
values n(Q), A(Q) or det(Q). Let L = conv{0, vi,..., vj} Q R d denote a <i-dimensional simplex and V 
a Delone triangulation. Then we consider the sets 

V L = {Q G 5> : |BRl(Q)| > 0}, 
V L = {QG^ :|BR L (Q)|=0}, 

w c = p|v L nA(z>). 

Lev 

Hence, by Corollary Oand Proposition|731the PQFs Q £ W v are those with /Jt(Q) < 1 in A(X>). Further, 
for D > and A > we consider 

Do = {Q G 5> : det(Q) > D}, 
M x = {Q £ Sy : X(Q) > A}. 

The following is known (see [DDRS70|, [Rys74]) and was partially discussed in previous sections: 

Proposition 9.1. Let L = conv{0, v\, . . . , v^} C M d be a d-dimensional simplex, T> a Delone triangula- 
tion, D > and A > 0. Then 

1. V^, and therefore Wp, is convex and bounded. 

2. Vl is a smooth (regular) hypersurface. 

3. Yin is unbounded, strictly convex and has a smooth boundary. 

4. Ma is unbounded and convex and has a piecewise linear boundary. 



5. Y)d and Ma are invariant with respect to the action of GL^Z) on S d 



>o- 



As a consequence of these properties, it is not hard to derive the following proposition (see [DDRS70|, 
|Rys74)). 

Proposition 9.2. Let V be a Delone triangulation. 
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1. The set of PQFs in A(D) attaining min Q(Q) is a single PQF, together with all of its positive 

Qe A(X>) 

multiples. It is invariant with respect to the action of D's automorphism group 

Aut(2?) = {U £ GL d (Z) : VU = V}. 



2. The set of PQFs in A(P) attaining min j(Q) is convex and contains a subset which is invariant 

QGA(D) 

with respect to Aut(P). 

A first, but non-geometric proof of property fl] was given by Barnes and Dickson [BD67]. There they 
also observed, as Ryshkov [Rys74] did for the lattice packing-covering problem, that it suffices to optimize 



among all PQFs whose automorphism group contains the group Aut(P). 

Corollary 9.3. If Q is a local optimum of the lattice covering or lattice packing-covering problem among 
all PQFs in A(D) whose automorphism group contains Aut(V), then Q is a locally optimal solution. 

Note that the statement above holds only if Q lies in the interior of the secondary cone. For local optima 
on the boundary of secondary cones we may apply the following trivial proposition. 

Proposition 9.4. A PQF Q is a locally optimal solution with respect to or 7 if and only if it is an optimal 
solution for all Delone triangulations V with Q £ A(£>). 

There probably exist two locally optimal solutions with respect to 7 in dimension 5 which both lie on 
the boundary of some secondary cones (see Section proof of Theorem II 1.2t . In dimension 5 and lower 
every locally optimal solution with respect to © lies in the interior of some secondary cone. However, in 
higher dimension there exist locally optimal solutions lying on the boundary. The Leech lattice gives such 
an example as recently shown by the authors in I SV04I . 



9.2 Rigorous Certificates 

Using convex optimization software to solve the covering or packing-covering problem, we are often limited 
to determining a certain certified range (see Section^ in which the optimum value is attained. To use Propo- 
sition |^]it is desirable to decide computationally whether or not the optimum is attained on the boundary 
or even on a specific facet of A(D). In many cases this is possible by using the following proposition which 
is a simple consequence of the convexity of determinant maximization problems. 

Proposition 9.5. Consider the determinant maximization problem 

minimize f(x) 

subject to G(x) y and F(x) t 0, K ' 

where f(x) = c l x — logdetG(a;). 

1. Let y 1 be the optimal solution of ©, and let y 2 be the optimal solution of © with the additional 
constraint F{x) h 0. If f( Vl ) < f(y 2 ), then F(y 2 ) = 0. 

2. Let y 1 be an optimal solution of © with the additional linear constraint n l x > 0, and let y 2 be an 
optimal solution of © with the additional linear constraint n l x < 0. 

(a) If/(y 1 )</(2/ 2 ) I thenn* 1/2 = 0. 

(b) If/( yi ) > /(y 2 ), then n*y 2 < 0. 
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Let Q pt denote a PQF with n(Q op t) = 1 attaining an optimum in A(£>) for B. The lattice packing- 
covering case is similar. Let I and u be a lower and an upper bound for a minimum of f(Q) = — log det(Q). 
Here and in the sequel we use nt>{Q) < 1 as an abbreviation for the LMI in Proposition 17.41 guaranteeing 
that the circumradius (with respect to Q) of all the simplices of the triangulation V is at most 1. 

The first part of Proposition 19.51 gives a sufficient criterion for Q op t lying on the boundary of A (£>). 
Suppose the upper bound of / on {Q G S d : Ht>(Q) < 1} is smaller than /. Then we have a certificate for 
Q opt G bd A(V). 

The first half of the second part of Proposition 19.51 gives a sufficient criterion for Q opt lying on a spe- 
cific facet F of A(D). Let denote a hyperplane containing F and let H~ , H + be the closed half- 
spaces containing A(D), respectively not containing A(D). If the upper bound of / on {Q G S d : 
Q G and ^v{Q) < 1} is smaller than /, then we have a certificate for Q op t £ i* 1 . 

The second half of the second part of Proposition 19.51 gives a sufficient criterion for Q opt not lying on 
a specific facet F of A(£>). If the lower bound of / on {Q G 5> : Q G and nv{Q) < 1} is larger 
than u, then we have a certificate for Q op t F. Clearly, such a certificate for all facets of A(D) gives a 
certificate for Qopt G A (2?). 

Note that these certificates give mathematical rigorous proofs when we proceed as described at the end 
of Section |6] 

9.3 Necessary and Sufficient Conditions 

Next we assume a PQF Q is given and we want to decide computationally if it is a locally optimal solution 
to the lattice covering or lattice packing-covering problem. In Section 19.3.21 we consider the case when Q 
lies in the interior of the secondary cone of some Delone triangulation. In this case we give necessary and 
sufficient conditions for Q being a locally optimal solution. The case when Q lies on the boundary of the 
secondary cone of some Delone triangulation is more subtle and we deal with it in Section l9.3.3l 

9.3.1 Normal Cones 

Before considering locally optimal solution we have to do some local analysis. Here, we shall compute the 
normal cones of the sets Wx>, and D£> at a given PQF Q. As a general reference to the basic concepts 
in convex and differential geometry used in the sequel we refer to the book [Sch93 1 by Schneider. 

We consider S d as Euclidean space with inner product (•, •). Let H = {S G S d : (N, S) = a} be a 
hyperplane with normal vector N G 5 a! \{0} and C a convex set with boundary point Q. Then H is said to be 
a supporting hyperplane of C at Q with outer normal vector N, if C C H~ = {S G S d : (N, S) < (N, Q}}. 
The normal cone of C at Q is then given by all outer normal vectors of supporting hyperplanes at Q together 
with the zero vector. Clearly the normal cones are convex. A hyperplane H is called a separating hyperplane 
for two convex sets with a common boundary point Q, if it is a supporting hyperplane of both sets, but with 
opposite outer normal vectors. Such a hyperplane exists if and only if the corresponding normal cones A/i 
and Af 2 at Q satisfy -Mi n N 2 ^ 0- 

Let us now compute the normal cones of the sets Wx>, Ma and Do at some PQF Q. 

Proposition 9.6. Let V be a Delone triangulation and Q G a PQF with fi(Q) = 1. Then the normal 

cone of Wp at Q is equal to — cone{gL(Q) ■ L G V and Q G Vl}, where ql{Q) = grad | BR^ \ (Q) 
denotes the gradient of the regular surface at Q. 

Proof. First, —gL{Q) / exists and is the unique outer normal vector of V l at Q, since is a regular 
surface defined by the polynomial equation |BRi(<5)| = 0. Second, we see that W©, in a sufficiently small 
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neighborhood of Q, is equal to the intersection C\g e ^ L \ L . Therefore, the normal cone of Wp at Q is 
equal to cone{— ql{Q) '■ Q £ ^l} (see ISch93l . Th. 2.2.1) and the assertion follows. □ 



By Proposition 19. II the boundary of D det (g\ is smooth at Q, hence up to scaling there is a uniquely 
determined outer normal. The following is well known (cf. Proposition 110.31 and its detailed proof in 
llW03l . Prop. 8.2.2). 

Proposition 9.7. The normal cone of D det (Q) at Q is given by — cone{Q -1 }. 

The normal cone of at a boundary point Q is determined by the shortest vectors of Q, that is, by 
those v € Ij d for which the homogeneous minimum A(Q) is attained. Recall that Q[v] = A for a fixed 
v G 7L d \ {0} is linear in the parameters qij, hence this condition defines a hyperplane in S d . 

Proposition 9.8. For v £ Z d \ {0} let V v G S d be the normal vector of the hyperplane {Q' G S d : Q'[v] = 
\{Q)} with (Q, V v ) = \{Q)- Then, the normal cone of Ma at Q is given by — cone{V v : Q[v] = A(Q)}. 

Proof. This follows immediately from 1Sch93l . Th. 2.2.1. □ 

Note, that a normal vector V v depends on the chosen inner product. For example for practical com- 
putations it is convenient to identify S d with R d ( d + 1 )/ 2 and to use the standard inner product (Q, V) = 

Y.i<j 1ij v ij and Vv = ( v ij) with v ij = ( 2 - 8ij)viVj. 

9.3.2 Interior Cases 

Because of the convexity of D det (Q), M^jqj and W-p a separating hyperplane at a PQF Q of with 
fi(Q) = 1 and one of the other two sets yields a necessary and sufficient condition for Q to be a locally 
optimal solution to either the lattice covering or the lattice packing-covering problem. Therefore a necessary 
and sufficient condition can be derived from the normal cones at Q with respect to these sets. 

Proposition 9.9. Let V be a Delone triangulation, and let Q £ A(2?) be a PQF with fi(Q) = 1. 

1. ( IBD67I ) Then, Q is a locally optimal solution to the lattice covering problem if and only if 

cone{Q -1 } n - cone{g L (Q) : L G V and Q G V L } ^ 0. 

2. Then, Q is a locally optimal solution to the lattice packing-covering problem if and only if 

cone {V v : (Q, V v ) = A(Q)} n - cone{g L (Q) : L G V and Q G V L } ^ 0. 

9.3.3 Boundary Cases 

Note that the foregoing proposition does not give criteria for PQFs Q on the boundary of some secondary 
cones. In such a case we need to replace — cone{<?£(Q) : L G V and Q G V^} by a generalized expression. 
That is done by considering for each Delone triangulation V with Q G A(£>) the normal cone Cx>(Q) 
of Wx) at Q. We have Ct>(Q) = — cone(A/i U A/2) with outer normal vectors Mi = cone{<?L(Q) : 
L G V and Q G V L }, as considered before, and outer normal vectors 

A/2 = {N G S d : (N, ■) = Q(l,L')(') f° r a P an " (L, L') of adjacent simplices of V with Q(l,L')(Q) = 0} 

from facets of A(D) containing Q. Since Dr> has a smooth boundary a necessary condition for a local 
covering optimum at Q is the convexity of Uqs a(P) ^ v at Q which is equivalent to PlQg A(i>) &>(Q) ^ 
0. In case of the lattice packing-covering problem, this convexity condition is not necessary for a local 
optimum. The following proposition however gives sufficient conditions for a local optimum in case the 
convexity condition is satisfied. 



23 



Proposition 9.10. Let Q G S d be a PQF with = 1. 

1. Then, Q is a locally optimal solution to the lattice covering problem if and only if 

cone{Q- 1 } n f] C V (Q) / 0. 

2. Then, Q is a locally optimal solution to the lattice packing-covering problem if 

cone{V v : (Q, V v ) = A(Q)} n f] C V (Q) + 0. 



Proof. The set Plgg ACP) ^>(Q) contains all common outer normal vectors of the sets Wp with Q G A(P) 
at Q. In particular, DQg A(X)) &>(Q) * s em pty if an d only if Uqe aTp) ^ p * s not convex at Q- ^ convex, 
then a separating hyperplane of Uq £ a(p) an< ^ ^det(Q)> respectively M A (q) at Q is sufficient for a local 
optimum at Q. The assertion then follows by Propositions I9.7l and 19.81 □ 



9.4 Conditions for Isolated Local Optima 

Locally optimal solutions to the lattice packing-covering problem may not be isolated optima, in contrast to 
solutions to the lattice covering problem. To determine computationally if a given locally optimal solution 
Q is isolated, we have to check if there exists a segment [Q, Q'] which lies in the boundary of M A (q\ and in 
the boundary of at Q. 

Proposition 9.11. Let V be a Delone triangulation and let Q G A(P) be a locally optimal solution to the 
lattice packing-covering problem with fi(Q) = 1. Then Q is not an isolated local optimum if and only if 
there exists an S G S d with 

1- (V v , S)>0 for all v G Z d with Q[v] = X(Q), 

2. (g L (Q),S) > for all L G V with Q G V L , 

3. h L {Q)S = for all L G V with Q G V L and (s L (Q), 5) = 0. 

Proof. The first condition says that there exists an e± > so that the segment [Q,Q + e\S] lies in Mwq) 
because M.\fQ) has a piecewise linear boundary. The first condition together with the second condition plus 
the local optimality of Q are equivalent to the fact that S is in the tangent space of at Q and in the 
tangent space of Mwq) at Q. Thus Q + e2<S is in the boundary of Wx> for sufficiently small €2 > if and 
only if the corresponding smooth hypersurfaces have curvature in direction S. This is equivalent to S 
being an eigenvector of eigenvalue of the Hessian hi,{Q) = hess | BR^ \(Q) which is the third condition. 
Hence, all three conditions are fulfilled if and only if there exists a segment [Q, Q'] in the boundary of the 
two sets M_wq) and W©, where Q' = Q + min(ei, 62) <5. □ 

Note again, that the condition applies only if Q is in the interior of a secondary cone. A similar sufficient 
condition for PQFs on the boundary of secondary cones can be attained by replacing ^^(Q) in 2 and 3 with 
the rays of the normal cone PlQg A(-p) Cx>{Q) of Uq £ a(P) ^ p at *3 ( see Proposition 19.101 1. if such exist. 

Finally, we propose two (as far as we know) still unanswered questions. 

Question 9.12. Does there exist a locally optimal solution to the lattice packing-covering problem which is 
not isolated? 
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Note that a positive answer to the following question would imply the existence of non-isolated locally 
optimal solutions. 

Question 9.13. Does there exist a (locally optimal) solution to the lattice-packing covering problem, which 
does not have all the symmetries of the corresponding Delone subdivision? 



10 Local Lower Bounds via Moments of Inertia 

In this section we give simple and efficiently computable local lower bounds for the lattice covering density 
and the lattice packing-covering constant. These bounds are "local" in the sense that they only apply to 
those PQFs lying in the topological closure of the secondary cone of a given Delone triangulation. For their 
computation we only need to know the coordinates of the simplices of the considered Delone triangulation. 
They are therefore useful tools in a heuristic search for "good" PQFs. The method goes back to Ryshkov and 
Delone. It is called the method of the moments of inertia because the central idea in its proof is analogous 
to the Parallel Axis Theorem of Steiner in classical mechanics. This analogy is explained in [ AM03 1. 

Let P C M. d be a finite set of points in (i-dimensional Euclidean space (M. d , (•,•)). We interpret the 
points of P as masses with unit weight. The moment of inertia of the points about a point x £ W 1 is defined 
as I X (P) = J2veP dist(ac 5 v) 2 . The centroid of P (center of gravity) is given by m = ^ ^2 v( =p v - From 
the equations 

dist(a3, i>) 2 = (x — v, x — v) 

= (x — m, x — m) + (m — v, m — v) + 2{x — m, m — v) 
= dist(a;, m) 2 + dist(m, v) 2 + 2{x — m, m — v), 

and X^ep( m ~ v ) = ®' we derive Apollonius' formula (see [Ber87| §9.7.6) 

I X (P) = \P\ distfo m) 2 + I m (P). (5) 

Hence, the moment of inertia of P about the centroid m is minimal. 

If the points of P form the vertices of a (i-dimensional simplex, then © gives a relationship between 
the radius of the circumsphere R, the center of the circumsphere c, and the moment of inertia about the 
centroid m of P: 



ir = v ' = dist (c,m) 2 + 



d+1 v ' d+1 

We can compute I m {P) using only the edge lengths of the simplex P. For every w € vert P we have 
by definition I W (P) = YlveP dist(io, v) 2 . Summing up and using © gives 



weP weP 



So, we get 



Im( - p ^ = dTT S dist (^,^) 2 - (6) 

{v,w}CP 



Let I? be a Delone triangulation of M. d , let L%, . . . , L n be the (i-dimensional simplices of the star of a 
lattice point (say for example the origin), and let rrij be the centroid of Li, i = 1, ... ,n. The arithmetical 
mean of the moments of inertia about the centroids of Lj with respect to a PQF Q is defined to be 



/7>(Q) = -X)WLi 
n — ■* 



n 

i=l 
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and is called the central moment of inertia of V with respect to Q. Note that we are now dealing with the 
inner product given by Q with dist(cc, y) 2 = Q[x — y] and that Iv{Q) is linear in the parameters qij of Q. 

Proposition 10.1. The central moment of inertia of V with respect to Q yields a lower bound for the 
inhomogeneous minimum of Q if V is a refinement of Del(Q). In this case we have 

KQ) > dTiW)- 

Proof. Let Ri be the radius and Cj be the center of the circumsphere of the simplex Lj, then 



H{Q) = max R 2 = max ( dist(cj, rrii) 2 + 

i=l,...,n i=l,...,n V 



d+1 

" * > • • • i ' 6 * 
— Cd+lln 

i=l 



□ 



To find lower bounds for the lattice covering density or for the lattice packing-covering constant, we 
have to minimize the linear function Ix> over all PQFs with a fixed determinant, or with a fixed homogeneous 
minimum, respectively. As in previous sections, we consider S d as a Euclidean space with inner product 

<v>- 

Proposition 10.2. There exists a PQF F G S^ with I v = (F, •). 

Proof. Since Id is a linear function, there is a F S 5 d with = (F, •). For every PQF Q, we have 

(F, Q) = I V (Q) > 0. Since S^ = {Q' G 5 d : (Q 7 , Q) > for all Q G 5^ }, meaning that is a 

'>o- 



self-dual cone, we have F G 5^ n . □ 



Because of the preceding proposition, minimizing the linear function Ix>, hence (F, •), over all PQFs 
with a fixed determinant D is geometrically equivalent to finding the unique PQF on the boundary of Drj — 
the determinant-D-surface — that has a supporting hyperplane with normal F (see Proposition I9.7I) . Thus 
we derive the following proposition (see [ Val03 1 Prop. 8.2.2 for a detailed proof). 

Proposition 10.3. A linear function /(•) = (F, •) with F G has a unique minimum on the determinant- 
D-surface. Its value is d\/ D detF and the minimum is attained at the PQF yf D detFF" 1 . 

Now we can plug Proposition 110.11 and Proposition 1 1 0. 3 1 to gether. This yields a local lower bound for 
the covering density of a PQF. 

Proposition 10.4. Let V be a Delone triangulation. Let Q be a PQF for which V is a refinement of Del(Q). 
Then we have a lower bound for the covering density of Q: 



d ^ d 



6(Q) > 6,(2?) = W [j^jj det F-K d , 

where F is the positive definite matrix given by the equation Iv(-) = (F •). Here and in the sequel we 
denote the local lower bound for the Delone triangulation V by 6*(D). 
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Proof. By Proposition 110.31 Id = (F,-) has the unique minimum d\/ D det F on the determinant-D- 
surface. Using this with D = det Q and applying Proposition 1 1 0. 1 1 we get the lower bound 6*(P) for 8(Q), 
because of 

(n(nu ,2 KQ) d ^ f Iv(Q) \ d lA d d detQdetF ( d \ d 

□ 

For a local lower bound on the lattice packing-covering constant we minimize the linear function 
Iv{-) = (F, •} over all PQFs with a fixed homogeneous minimum. In analogy to the above, we replace 
the determinant- £>-surface (the boundary of Dp) by the homogeneous-minimum-A-surface (the boundary 
ofM A ). 

Proposition 10.5. Let A be a positive number. Let V be a Delone triangulation and let Q be a PQF for 
which V is a refinement of Del(Q). Further, let the minimum of the linear function Id(-) = (F, •) on 
the homogeneous-minimum-A-surface be attained at a PQF Qf,x- Then we have a lower bound for the 
packing-covering constant of Q: 



7 (Q) > 24 



(F,Q 



F,X) 



(d + l)A 



Proof. The PQF ^jp-QF,\ attains the minimum of (F, •) among PQFs in Muq-). Applying Proposition llO.il 
gives the desired lower bound: 

(7(Q)/2) " HQ) - {ITT) ,KQ) ~ (d+l)-A(Q) - TdTI)A • 

□ 



It is difficult in practice to obtain Qf,x and therefore the lower bound of the proposition. Instead, if we 

v 

X 



minimize I-p on an approximation of M.\, suitable for V, things become easier. We set 



QeS d : Q[v] > A for all v e Z d with [0, v] G p| . 



Now finding the minimum of I-p on M^ 1 reduces to a linear program with finitely many constraints. By 
looking at equation © it is easy to see that this linear program is bounded from below since every summand 
is at least A. Because of M\ C we obtain the following practically useful lower bound for 7. 

Proposition 10.6. Let A be a positive number. Let V be a Delone triangulation and let Q be a PQF for 
which V is a refinement of Del(Q). Further, let the minimum of the linear function I-d(-) = (F, •) on 
be attained at a PQF Q?,. Then we have a lower bound for the packing-covering constant of Q: 



i(Q)>i*{v) 



' (FQ3,x) 
(d+l)-X' 
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11 Computational Results 



In Section[8]we developed algorithms for the solution of the lattice covering and the lattice packing-covering 
problem. Here, we want to demonstrate that these algorithms are not purely of theoretical interest. We 
implemented the algorithms in C++, using the package MAXDET 3 of Wu, Vandenberghe and Boyd and the 
package Irs 4 of Avis as subroutines. The interested reader can download the implementation from our web 
page MM- 

With the implemented algorithm we are able to determine the solution of the lattice covering problem 
and the lattice packing-covering problem in dimensions d = 1, . . . , 5. This hereby reproduced the known 
results. For dimension d = 5 this computation took on a 2GHz Intel Pentium computer less than 90 
minutes. By our computations we get a conjectural list of approximations of all locally optimal solutions. 
In Section riT.il we give some details. 

More important, with the help of a heuristic method, which we describe in Section [T 1.21 this approach 
produces interesting new lattices in dimensions d = 6 which give better coverings and packing-coverings 
than previously known ones. We analyze them in Section [T 1.31 and in Section IT 1 .41 where we give rigorous 
proofs of some facts concerning these lattices which we found experimentally. This analysis turned out to 
be fruitful: The new 6-dimensional lattices were our starting point to find new best-known covering lattices 
in dimensions 7 and 8. We describe the 7-dimensional case in Section II 1.51 The 8-dimensional case is 
explained in ISV04I . 

11.1 Dimensions 1, . . . , 5 

We described the solutions to both problems for the 1, 2, and 3-dimensional cases in Section |3] and |4] Since 
in these cases there is only one type of Delone triangulation, these results are computationally rather trivial. 

The 4-dimensional case requires more work because there are three inequivalent Delone triangulations 
in dimension 4. Baranovskii [Bar65 1 and independently Dickson [Dic67 1 found all three inequivalent locally 
optimal solutions to the lattice covering problem in dimension 4. Earlier Delone and Ryshkov showed that 
the lattice A| gives the best 4-dimensional lattice covering without determining all locally optimal solutions. 
We describe their approach in Section [Tl.2l 

For the lattice packing-covering problem the 4-dimensional case was resolved by Horvath [Hor82|. He 
shows that there exist, as in the covering case, three isolated, locally optimal solutions — one for each 
Delone triangulation. It is interesting that Horvath's optimal packing-covering lattice H04 does not belong 
to the family of root lattices and their duals. An associated PQF is 
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The first summand, associated to the best packing lattice D4, lies on an extreme ray of the secondary cone 
belonging to Q Ho4 . 

In a series of papers Ryshkov and Baranovskii solved the 5-dimensional lattice covering problem. In 
|Rys7 3 1 Ryshkov introduced the concept of C-types. Two Delone triangulations are of the same C-type 
if their 1-skeletons (the graph consisting of vertices and edges of the triangulation) coincide. He gave an 
algorithm to find all inequivalent C-types in any given dimension. He computed that there are 3 inequivalent 
C-types in dimension 4 and that there are 76 inequivalent C-types in dimension 5. Using this list Baranovskii 

; http: //www. Stanford . edu/ ~boyd/MAXDET . html| 
' http : / / cgm. cs .mcgill ■ ca/ ~ avis /C/ Irs . html| 
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and Ryshkov enumerated 221 (of 222) inequivalent 5-dimensional Delone triangulations in [BR73]. They 
described the triangulations in more detail in [BR75|. In the last paper of the series [RB75| they showed 
that the lattice A5 provides the least dense 5-dimensional lattice covering. In their proof they do not find 
all locally optimal lattice coverings. By using estimations (the method explained in Section [T 1.21 is one of 
their main tools) they merely show that all local minima exceed the covering density of A5. As mentioned 
in Section 1531 Ryshkov and Baranovskii missed one Delone triangulation. Fortunately, this does not give a 
thinner lattice covering than A5. 

Using our algorithm, and the techniques described in Section to obtain certificates we produced a list 
of inequivalent locally optimal lattice coverings in dimension 5. This list is conjecturally complete and gives 
approximations of the local optima. The computation, including all certificates, takes about 90 minutes on a 
2GHz Intel Pentium computer. We can prove the following theorem rigorously. 

Theorem 11.1. In dimension 5, there exist at least 216 and at most 218 inequivalent, local minima of the 
lattice covering density function ©, ranging from approximately 2.124286 to approximately 2.757993. All 
of them are attained in the interior of their secondary cones. 

Proof. We take a complete list of 222 inequivalent secondary cones, generated by sec (see Section [5] and 
our web page [SV05|). For each cone, we applied the program coop {covering optimizer), which is based 
on rmd {rigorous MAXDET) and available on our web page [SV05|. It computes certified bounds for © 
within the closure of a given secondary cone by solving the determinant maximization problem described in 
Section l8~T1 The bounds for the duality gap described in Section[6]are mathematical rigorous, since we use 
rational arithmetic only. Solving additional determinant maximization problems as described in Section l9~2l 
the program then tries to obtain certificates for the approximated local optimum to be attained in the interior 
or on the boundary of the secondary cone. As a result, there are 216 certified local optima attained in the 
interior of secondary cones and 4 certified non-optima attained on the boundary. This leaves two cases 
which might be local optima, attained in the interior of their secondary cones. However, if they are attained 
on the boundary (and the numerical evidence strongly supports this) then they are not local optima, since 
the corresponding secondary cones adjacent to the facets in question have a smaller local minimum. □ 

Note that our computations prove rigorously that the PQF of Barnes and Trenerry [BT72] yields the 
second best locally optimal lattice covering with density of approximately 2.230117. 

In 1986, Horvath [Hor86] solved the lattice packing-covering in dimension 5. The proof is about 70 
pages long (private communication). As in dimension 4, Horvath's lattice H05 is not among previously 
known ones. An associated PQF is 

-2 -2 -2\ 
-2 -3 -2 
6 -1 -1 . 
-16 
-10 6 / 

Again, the first summand is associated to the best packing lattice D4. It has rank 4 and lies on an extreme 
ray of the secondary cone belonging to Qho 5 - 

By applying propositions of Section[9]we were able to reproduce Horvath's result and moreover to attain 
a list of locally optimal solutions. Again this list is conjecturally complete and gives approximations of the 
local optima. Compared to the covering problem, there are many more secondary cones that contain no 
locally optimal solution. The computation, including all certificates, also takes about 90 minutes on a 2GHz 
Intel Pentium computer. We can prove the following theorem rigorously. 
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Theorem 11.2. In dimension 5, there exist at least 47 and at most 75 local minima of the lattice packing- 
covering constant 7, ranging from approximately 1.449456 to approximately 1.557564. At least 45 of them 
are attained in the interior of their secondary cones. 

Proof. As for the proof of Theorem lll.il we used the complete list of 222 inequivalent secondary cones, 
generated by sec. For each cone, we applied the program pacoop (packing-covering optimizer), which 
is also based on rmd and available on our web page HSV05I . It computes certified bounds for 7 within the 
closure of a given secondary cone by rigorously solving the semidefinite program described in Section IOI 
Again the program tries to obtain certificates for the approximated local optimum to be attained in the 
interior or on the boundary of the secondary cone. As a result, there are 45 certified local optima attained 
in the interior of secondary cones and 147 certified non-optima attained on the boundary. There are 30 of 
222 remaining cases in which the program pacoop did not give a certificate for the optimum to be attained 
in the interior or the boundary of the secondary cone. However, if all of them are attained on the boundary 
(and the numerical evidence strongly supports this), there are two local optima on the common boundary of 
three inequivalent secondary cones each. Thus there exist at least 47 local optima as claimed. □ 

11.2 Heuristic Methods 

Before going to dimension 6, let us explain one important heuristic method, which is essential for finding 
the new lattices. One of the biggest problem in finding good lattice coverings or packing-coverings in 
dimension 6 and higher is that it is not apriori clear which Delone triangulations admit good ones. Solving a 
determinant maximization problem is a rather time consuming task and there are a lot of inequivalent Delone 
triangulations to consider. So a desirable tool is a fast computable lower bound, as described in Section fTOl 
We view the set of Delone triangulations as an undirected labeled graph. A node represents a Delone 
triangulation and two nodes are adjacent if their Delone triangulations are bistellar neighbors. Let V be a 
Delone triangulation. We label its node by the local lower bound ©*(£>) or 7* (£>). We can use the labeling 
in two different ways. 

On the one hand it is clear that if the labeling of a node is large, then the considered Delone triangulation 
does not admit a good lattice covering and lattice packing-covering respectively. Ryshkov and Delone 
[DR63| solved the lattice covering problem in dimension 4 by this method. However, the five-dimensional 
lattice covering problem cannot be solved in this way. From the 222 inequivalent Delone triangulations 
there are 20 whose local lower bound is smaller than ©(A^). 

On the other hand we can hope that V admits a good lattice covering or packing-covering if the local 
lower bound is small. In dimension 6 the hope that good local lower bounds yield good lattice coverings is 
partially fulfilled. We report on a typical example: We started from the Delone triangulation of Voronoi's 
principal form of the first type and took a random walk of length 50. Then, we found a node labeled 
by ©*(D) ~ 2.585149. From this we proceeded by taking a neighboring node having the smallest local 
lower bound. By repeating this greedy strategy we resulted in a node labeled by ©*(£>) ~ 2.318034. This 
node is interesting for several conjectural "extremeness" properties. It yields the smallest known local lower 
bound and it has the largest known number of neighbors, namely 130. As we will see in the next section there 
exists a locally optimal lattice covering which belongs to this node with covering density of approximately 
2.466125. At present this is the second best known 6-dimensional lattice covering which is locally optimal. 
Furthermore we will see that there exists a locally optimal lattice packing-covering belonging to this node 
which currently defines the best known 6-dimensional lattice packing-covering. 
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11.3 New Six-Dimensional Lattices 



We have not been able to solve the lattice covering problem in this dimension. However we found some 
new interesting covering lattices. As we reported in Section |3j Ryshkov [Rys67] asked what is the first 
dimension d where does not give the least dense lattice covering. 

A classification of all inequivalent Delone triangulations in dimension 6 is not within reach. So we can 
not use our algorithm to find the best lattice covering. Nevertheless we ran it partially to find good lattice 
coverings. Since there are no other good lattice coverings in the neighborhood of Voronoi's principal form 
of the first type, we used the heuristic method described in the previous section. In this manner we found 
about 100 inequivalent secondary cones containing lattice coverings better than Ag. Thus, we can give an 
answer to Ryshkov's question: 

Theorem 11.3. Dimension d = 6 is the smallest dimension where the lattice does not give the least 
dense lattice covering. 

Proof. This theorem follows from Theorem lll.7l or from Theorem ll 1.81 □ 

Two of the new lattice coverings were strikingly good. Of course, by the computational optimization we 
only got numerical estimates of these coverings. 

The now second best-known lattice covering which is locally optimal was quite easy to find: running 
our heuristic method described above finds this lattice covering in most of the trials. An approximation is 
given by the PQF 
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and its covering density is @(Qq 2 ) 2.466125. The Delone subdivision is a triangulation and its secondary 
cone has 130 facets. The local lower bound is approximately 2.318034. 

For a while we thought that this might by the best lattice covering in dimension 6, but then 
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with covering density Q(Qq 1 ) rj 2.464802 came up. The Delone subdivision is a triangulation and its 
secondary cone has 100 facets. The local lower bound is approximately 2.322204. After this, we did not 
find any further lattice covering records in dimension 6. Furthermore, we did not find a 6-dimensional 
Delone triangulation whose secondary cone has more than 130 facets or whose local lower bound is less 
than 2.318034. 
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Nevertheless, in the secondary cone of Qf we found the PQF 
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which gives currently the best known lattice packing-covering in dimension 6 with packing-covering con- 
stant 7(Qg C ) ~ 1.411081. In the next section we will examine these new PQFs in greater detail. 



11.4 Beautification and a Unified View 

Although we found an answer to Ryshkov's question, these results are not fully satisfying. We want to 
know the exact lattices and prove rigorously that they have a good covering density and that they are locally 
optimal. Even more important, we want to know an interpretation of why these lattice coverings are good. 
To accomplish this we collect some more data. 

The automorphism group of Del(Qg 2 ) has order 3840 and the one of Del^g 1 ) has order 240. With 
the knowledge of the groups we were able to compute the extreme rays of both secondary cones [DV04|. 
The secondary cone Ci = A(Del(Qg 2 )) has 7, 145, 429 and the cone C 2 = A(Del(Q§ 1 )) has 2, 257, 616 
extreme rays. Both contain an extreme ray associated to the lattice Eg given for example by the PQF 
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After transforming Qq and Qq by integral unimodular transformations we can assume that the Delone tri- 
angulations of the two PQFs are refinements of the Delone subdivision Del (QE^)- This Delone subdivision 
was investigated in different contexts (see [Wor87|, [CS91|, [MP95| and [Bar92]). We briefly review the 
main results: 

Proposition 11.4. In the star of the origin are 720 full-dimensional 6-dimensional Delone polytopes, and 
the automorphism group of Qe* acts transitively on these full-dimensional polytopes. Each polytope is the 
convex hull of three regular triangles lying in three pairwise orthogonal affine planes. Each polytope has 9 
vertices, 27 facets, and three different triangulations, where each triangulation consists of nine 6-dimensional 
simplices. The covering density is Q(Qe*) = • ^6 ~ 2.652071. 

Proof. Except for the possible refining triangulations, all this data is well-known, see for example [CS91 1, 
Summary for Eg. 

To describe the triangulations we introduce coordinates. A Delone polytope of Qe* is similar to the 
polytope P = convlu!,^,^,^!,^,^,^!,^,^, } where u x = e x , u u = -\e x + ^e 2 , Uu = 
-\e\ - ^e 2 , f i = e 3 , v w = -±e 3 + ^e 4 , Vjj = -±e 3 - ^e 4 , w x = e 5 , = -|e 5 + ^e 6 , 
Wjj = — — -^eg. The three possible triangulations are given by the set of nine 6-dimensional simplices 

T u = {convvert P \ {v z , w z >} : z, z G {1, u;,o;}}, 
T v = jconv vert P \ {u z , w z /} : z, z' £ {1, u>, a;}} , 
T w = {conv vert P \ {u z , v z '} : z, z' G {1, uj, a;}} . 
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To finish the proof one has to show that these sets indeed define triangulations and that they are the only 
possible triangulations. This can be done by a straightforward computation using the facts that the minimal 
affine dependent subsets of vert P are 

{v 1 ,V^,Vjj,W 1 ,W^,W I j},{u 1 ,U uj ,U 7 j,W 1 ,W uj ,W Z j},{u 1 ,U uj ,U Z j,V 1 ,V u) ,V 7 j}, 

that the facets of P are convvertP \ {u z ,v z i,w z »} where z,z',z" G {1,cj,uJ}, and by applying the 
following proposition. 

Proposition 11.5. ([Ram97 1, Prop. 2.2) Let A C R d be a finite point set, and let T be a set of d-dimensional 
simplices with vertices in A. The set T defines a triangulation of the polytope conv A if and only if the 
following two conditions hold: 

1. For all S,S" 6 T there exists a minimal affine dependency J2 a eA A a a = with J2 a eA A a = so 
that {a G A : A a > 0} C S and {a G A : X a < 0} C S'. 

2. For all S G T and for every (d— 1) -dimensional facet F of S there exists either a (d — 1) -dimensional 
facet F' of conv A with F C F' or there exists another simplex 5' G T \ {5} also having F as a 
facet. 

□ 

Using this information we are able to prove that the PQF Qq 1 is closely related to Qe* ■ 

Theorem 11.6. The PQF Qq 1 gives the least dense lattice covering among all PQFs whose Delone subdivi- 
sion is a refinement of the Delone subdivision Del(QE*)- 

Proof. Our proof is computational and uses a branch and cut method. We have to show that all secondary 
cones of a Delone triangulation refining Del(QE*) do not contain a PQF with covering density less than 
G(Qg 1 ). There are 40 full-dimensional Delone polytopes of Qe* which cannot be transformed into each 
other by translations or by the map x h-> —x. Since each of these Delone polytopes has 3 possible trian- 
gulations, the number of all periodic triangulations refining Del(QE*) is 3 40 . It is not apriori clear how to 
distinguish between Delone and non-Delone triangulations, and it is not possible to generate all 3 40 triangu- 
lations. We choose a backtracking approach instead. 

We arrange partial triangulations in a tree. On every level one of the 40 Delone polytopes is triangulated 
so that we have 3 n nodes on the n-th level. For every node N we define the value 

@N = max |det(<5) : Q G 5> , BRl(Q) >z for all simplices L of partial triangulation iV} . 

Obviously, G^r is a lower bound for the covering density of any PQF whose Delone subdivision refines the 
partial triangulation N. We can compute a lower bound of this value by solving a determinant maximization 
problem similar to the one in Section 18.11 Note that this can be done rigorously using rational arithmetic 
only if we proceed as described in Section[6] If the lower bound is larger than 2.464802, we can cut the tree 
at this node, since the covering density of Qq 1 is less than 2.464801 (see Theorem ll 1.81 . 

This algorithm visits exactly 432 nodes of depth 40. One of these triangulations equals De^Qg 1 ). All 
the others are equivalent to De^Qg 1 ) because the automorphism group of De^Qg 1 ), which has order 240, 
is a subgroup of Aut(Q E * ) and the order of Aut(Q E * ) equals 103680 = 432 • 240. □ 

This computational proof takes about two weeks on a 2GHz Intel Pentium computer. The source code 
e 6d . cc is available from our web page [SV05 1 as part of the package rmd. 
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The knowledge of the automorphism groups of De^Qg 1 ) and Del(Qg 2 ) also enables us to give a unified 
view on both lattices. We have 

Aut(Del(Qg 1 )) C Aut(Del(Qf )) C Aut(Q E *)- 



The automorphism group Aut(Del((5§ 2 )) turns out to be the subgroup of Aut(QE*) stabilizing the mini- 
mal vectors ±ei, and Aut^Del^g 1 )) is the intersection of the two subgroups of Aut(QE*) stabilizing the 
minimal vectors dbej and ±e2 respectively. 

The subspace Ii of all quadratic forms invariant under the group Aut(Del(Qg 2 )) is spanned by the 
PQFs Qe* and Ri (see below). At the same time, Ii n Ci is a cone with extreme rays Qe*,Ri- By 
Proposition Qg 2 has to lie in cone{QE;t ,R\}- The subspace 1 2 of all quadratic forms invariant under the 
group Aut(Del(Qg 1 )) is four-dimensional. The cone I2 D C2 has six extreme rays Qe*,R2, ■ ■ ■ , R&, where 
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Note that R2 lies in \\. Altogether, this yields the "picture" in dimension 21 given in Figure 5. 

Let us finally try to find the exact coordinates of the PQFs. This is easy for Qq 2 . We know that we can 
scale Qq 2 so that Qq 2 = Qe* + xR\, for some x £ M>o- Now, the exact finding of x boils down to finding 
roots of an univariate polynomial: we have to minimize the function x 1— ► /u(Qe* +xR\) d / det(QE* +xR\) 
where we know, because of the approximate solution, that h(Qe* + xR\) is a polynomial for all points in a 
sufficiently small neighborhood of the exact x. This leads to the Ansatz 

V1057 - 1 
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Now we can use the tools we introduced in Section[9]to prove rigorously that this Ansatz works. 
Theorem 11.7. The PQF Qq 2 gives a locally optimal lattice covering. Its covering density is 
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Kg 2.466121650. 



Proof. Our proof is again computational. We provide the MAGMA program check_q6c2 . m at the 
arXiv.org e-print archive. To access it, download the source files for the paper math.MG/0403272. 
Since it uses rational arithmetic only, the proof is rigorous. 

L •■ \ r'l 1 I-.. Ill II I ' WIM / ». I" / ^ 



Let us describe the steps. First we compute the Delone subdivision of Qq 2 : We compute the Delone 



il = To^Ee + TO^l' ^2 = tqVE* f Tom, 
Then we check that they coincide and that it turns out to be a triangulation. 



subdivisions of the three PQFs Q 1 = j^Qe* + jqRi, Q2 = j^Qe; + jqRi, and Q 3 = ±lQ E ; + jqRi- 



Then we show that Qq G 
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Figure 5. Unified view on Q E *, Qf,Qf and Q^ c . 



convjQi, Q2, Q3} so that Del(Qg 2 ) = Del(Qi). Now we can compute the circumradii of all the Delone 
simplices by the formula © in proof of Proposition l7.il This gives the value 0(Qg 2 ). Finally we prove that 
(Jq 2 gives a locally optimal lattice covering using the criterion of Proposition |^J We compute the gradients 
gi for all simplices L 6 Del((5g 2 ) achieving maximum circumradius. Summing them up yields a multiple 
of-iQf)- 1 . " □ 

This computational proof takes about one minute on a 2GHz Intel Pentium computer. 

What is the general pattern behind the beautincation process? Let Q be a locally optimal lattice covering 
with Delone triangulation V. We use the symmetry of V to find the subspace in which Q lies. This reduces 
the number of involved variables. The simplices of the Delone triangulation which have circumradius 1 give 
equality constraints. Then we maximize the determinant of the quadratic forms lying in the subspace subject 
to the equality constraints. For this optimization problem, which involves only algebraic equations, we can 
use Grobner basis techniques. 

Unfortunately, we were not able to solve the corresponding algebraic equations for Q^ 1 so we have to 
be satisfied with an approximation. 

Theorem 11.8. The covering density of the PQF Qf is at most 2.464801. 

Proof. The covering density of the PQF Q E * + xR 2 + yR 3 + zR A with x = 0.15266328480099, y = 
0.32884740614948, z = 0.13827491241153 is smaller than 2.464801. For a computational proof of 
this fact we provide the MAGMA program check_q6cl . m, available from the source files of the pa- 
per math .MG/0403272 at the arXiv . org e-print archive. The steps are similar to the first steps of the 
proof of Theorem ll 1.71 □ 

On the basis of Theorem 111.61 and our extensive computational experiments we make the following 
conjecture. 

Conjecture 11.9. The PQF Qq 1 provides the unique least dense lattice covering in dimension 6. 
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Using a similar and more successful beautification process for Q^ c , we make the Ansatz 



v'798- 18 
Q e - Qe* + ^ Hi. 



Theorem 11.10. Qq C is a locally optimal solution to the lattice packing-covering problem, lying in the 
interior of its secondary cone. Its lattice packing covering constant is 



7 (QP C ) = 2y 2V798 - 56 ss 1.411081242. 

Proof. This is similar to the proof of Theorem II 1.71 We provide the MAGMA program check_q6pc . m, 
available from the source files of the paper math.MG/0403272 at the arXiv.org e-print archive. □ 



11.5 Dimension 7 

After analyzing the 6-dimensional case, we got a feeling of where we have to search for good 7-dimensional 
lattice coverings. We took Qe* and a lattice vector which is a longest vector of the shortest vectors in the 
cosets Z, d /2Z d . We computed the stabilizer group of this vector and the invariant subspace of this group. 
By perturbing Qe* in this subspace randomly, we found a PQF whose Delone subdivision is a triangulation. 
We solved the determinant maximization problem of Section l8~Tl which belongs to this Delone triangulation 
and found the PQF 
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We are quite surprised that this PQF has rational entries. 

Theorem 11.11. Qj is a locally optimal solution to the lattice covering problem, lying in the interior of 
its secondary cone. Its inhomogeneous minimum is \i = its determinant is det Qy = 2 • ll 6 , so that 
Q(Qy) ps 2.900024. 

Proof. Again, this is similar to the proof of Theorem ll 1.71 We provide the MAGMA program check_q7c . m, 
available from the source files of the paper math.MG/0403272 at the arXiv.org e-print archive. □ 
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